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Abstract 

A new method for predicting chemical rate constants using thermodynamics has been applied 
to the hydrogen/oxygen system. This method is based on using the gradient of the Gibbs free 
energy and a single proportionality constant D to determine the kinetic rate constants. Using 
this method the rate constants for any gas phase reaction can be computed from 
thermodynamic properties. A modified reaction set for the H/O system is determined. All of 
the third body efficiencies M are taken to be unity. Good agreement was obtained between 
the thermodynamic method and the experimental shock tube data. In addition, the hydrogen 
bromide experimental data presented in previous work (ref. 1) is recomputed with M’s of 

unity. 

Introduction 

There are several reasons why it is desirable to relate the reaction kinetic rate constants to 
thermodynamic data. First, modem kinetic models contain many reactions with many 
empirical constants to form large mechanisms. These systems are not universal and are 
complex to assemble from the experimental data base. In contrast, thermodynamic equilibrium 
calculations have been simplified and systematized and are only species dependent. If the 
kinetic constants can be related to thermodynamics, it would significantly reduce the work 
involved in doing kinetic calculations. 

Secondly many kinetic constants are still unknown; if these reactions could be computed by 
using thermodynamic data, the kinetics calculation could be done without performing more 
experiments. 

There already exist several relationships that link thermodynamics to the kinetic rates. In this 
report an expression is derived which links the individual rate constant to the thermodynamic 
data. For instance in kinetics textbooks: the ratio equation (1) relates thermodynamic 
equilibrium constant K,. to the ratio of the forward rate constant lq to the reverse rate constant 
k,: 
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( 1 ) 



Normally, one thinks of kinetics as a path process, however the kinetic rate constants are state 
variables like K,.. The kinetic constants k are functions only of temperature and are therefore 
state variables and not path variables. Equation (1) holds not only at equilibrium, but also at 
every point of the reaction from initial conditions to final state. Furthermore, it is 
independent of concentration. 

It is widely thought that an overall relationship between kinetics and thermodynamics is not 
possible. All of the kinetic text books principally state that thermodynamics is only for 
equilibrium calculations whereas kinetics describes the rate of approach to equilibrium; thus, 
they state that thermodynamics would not be useful for kinetics calculations. However, as 
stated, relationship (1) holds for all time, not only at equilibrium. 

With thermodynamics, the complete energy surface is available for ideal gas mixtures for all 
compositions (from initial reactants to final products), not only for equilibrium states. This 
information can be used to compute the kinetic rate constants. The change of system free 
energy with reaction dG/dx. is known, and is used to compute the kinetic rate of reaction. The 
current use of system free energy is that when the free energy is minimum, the system is in 
equilibrium and no further reaction can occur. When the system free energy is above the 
minimum, reaction can occur. The Gibbs free energy is definitely the driving force for the 
reaction. We state that the reaction rate is related the the free energy gradient using a 
proportionality constant D between free energy and the chemical rate. Then using this 
relationship, the individual kinetic rate constants can be calculated. Only one number for D is 
used for the complete regime of temperature, pressure, and composition. 

The following discussion compares the classical method with the thermodynamic method for 
the H/O system. In doing this a third body efficiency M needed to be included in the new 
method to obtain the correct pressure behavior. M is set to unity for all molecules. A 
comparison of the two methods is made with experimental data for the H/O system. Finally 
the HBr results from reference 1 are recomputed including M taken as unity and making all 
stoichiometric coefficients into integers rather than using fractions. In general the agreement 
with data is very good. 

The Classical Method of Kinetic Rate Constants for the H/O System 

The classical kinetic method can be described for a particular rate by: 

V 1 S 1 +V 2 S 2 = V 3 S 3 +V 4 S 4 > (2) 

where the v. are the stoichiometric coefficients for species Sj. The forward 
rate r f is given by: 
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(3) 


r f - kf Ci C 2 . 

and the kinetic rate constant kf is defined as: 

k f = Af T n exp(-E/RT) (4) 

The kinetic constants A f , n, and E are determined by the best fit of the experimental data. A 
large data base is used to fit the many rate constants. 

The data of reference 2 is used for comparing both the clasical and thermodynamic methods. 
This data is also published in references 3 and 4. This data follows the OH concentration 
using laser absorption in a shock tube experiment. Weak mixtures of H 2 and 0 2 in Argon are 
used to try to isolate the principle reactions. The range of test conditions was wide from 
equivalence ratios of 0.05 to 10., pressures of 0.7 to 4.1 atm, and temperatures from 1050 K 
to 2500 K. Over 200 test points are reported. 

Theclassical mechanism given in Reference 2 is shown in table 1 : 


Table 1. Classical Rate Equations from Ref. 2. 
Name Reaction 

CL1 H +02 = OH + O 
CL2 O +H2 = OH + H 
CL3 OH +H2=H20+H 
CL4 O + H2O=2.0OH 

CL5 0 + 0 = 02+M 

AR=1., H2=2.9, 02 = 1.2, H20=18.5 
CL6 H + H = H2+M 

AR=1., H2=4.0, H=26., H20=12. 

CL7 H + O = OH +M 
AR=1., H2O=5.0 
CL8 H + OH = H20 + M 

AR=1., H2=2.5, H20= 16.25 
CL9 H + 02 = H02 + M 

AR= 1., H2=3.33, 02 = 1.33, H20=21.3 
CLIO H02+ H = OH +OH 
CL11 H02+ H = H2 +02 
CL12 H02+ H = H20 + O 
CL13 H02+ O = 02 +OH 
CL14 H02+ OH =H20+02 
CL15 H02+ H02=H202+02 
CL16 H2O2+M=2.0 OH+M 

AR=0.67, 02=0.78, H2O=6.0 
CL17 H202+H=H02+H2 
CL18 H202+H=H20+0H 
CL19 H202+0=H02+OH 
CL20 H202 + OH = H20 + H02 



3 



The classical reaction mechanism consists of 20 reaction steps. Each reaction has three 
experimental constants with a total of 14 third body efficiencies for a total of 74 constants. In 
this system there are 9 species, H 2 , H, 0 2 , O, OH, H z O, HO z , H 2 0 2 and Argon. Reactions 
CL1, CL2, and CL3 are the principle reactions with CL1 and CL2 producing OH and then OH 
reacting with H 2 through CL3 to produce H z O. 

The differential kinetics equations were integrated with the NASA Lewis chemical kinetics 
code LSENS, reference 5. The path was taken to be a constant density case to represent the 
shock tube data of reference 2. LSENS was used for both the classical and thermodynamics 
methods. The thermodynamic properties for both the classical and thermodynamic methods 
were obtained from reference 6. The thermodynamic entropy reference b2 was corrected from 
one bar pressure to one atmosphere by: ft2(lof/w)=ft2(lAar)-0.01316 for all species. 

Thermodynamic Method of Computing the Kinetic Rate Constants 

In the thermodynamic derivation, the kinetic rate is assumed to be related to the gradient of the 
free energy with respect to the reaction progress variable x, or: 

rt “ e><P RT ftT (5) 


The free energy G , ref. 7, is given by: 


G=E ("i RT ln(p.)) 

i i 


Differentiating equation (6) with respect to x, yields: 


1 dG_y^ 
RT 1h j 


RT 


^-E Vri 


o 

i 


RT 


+ E( V fi ln Pfi)-E( V ri to Pri> 


( 6 ) 


( 7 ) 


The standard thermodynamic expression for Kp is the first two terms on the right hand side 
and is expressed as: 


K p = exp 


v fi Pfi 
RT 


-E 


0 

i 


RT 


, -AG°* 
or = exp( ) 

RT 


( 8 ) 


Kp is independent of composition or degree of reaction. It represents the difference in 
chemical potential between the reactants and products. In addition K,. is given by (where 
R=82.05 atm-cc/mole-K): 


K c =K p (RT) AVi =K p (RT) (2Vf ‘ 2w ' ^ . 


( 9 ) 


At this point, it is important to distinguish chemical kinetic constants, kp and k,.. When the 
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kinetic rate is given in terms of concentrations, k,. is used; and when the rate is given in terms 
of pressure p, kp is utilized. 


cX? ; 4>=k p,V ; and C,=A 

c 1 2 > Q rl rj, y t TiT 1 


Then: 


dt 

k e =k p (RT) 1 


RT 


,(Sv ti -D . 


( 10 ) 


similar to K,. and Kp in equation (9). 

Returning to equations (5) and (7) and equating the pressure terms (IE. the 
concentration terms) to the concentration terms in equation (3) for the forward rate, one 
obtains: 


k cf =DK p (RT) SVfi l for K p < 1. (11) 

where D is a constant to be determined. Equation (11) provides the connection between 
thermodynamics and the forward chemical kinetic constant that shall be used. 

Equation (1 1) without the Kp inequality cannot hold for both the forward and reverse rates 
because equation (1) would not be satisfied. Remember that Kp^l/Kpf, so when Kp f < 1, then 
K pr > 1. To satisfy the kinetic constant ratio equation (1), the reverse rate must be given by: 

k =D(RT) Sv ”' 1 for K > 1. (12) 

CT pi 


The proportionality constant D can be a universal constant. But as will be shown in this report, 
a different value had to be used for the H/O system, D=7.0xl0 7 , compared to the HBr 
computation where D=3.0xl0*. For all of the H/O calculations, D was held constant at 7xl0 7 
for all the reactions. The difference in the values for the two systems could be the result of 
impurities in the HBr measurements which were made in 1906. No sensitivity to surface area, 
air addition or water vapor was found; but the system was sensitive to iodine which produced 
pronounced inhibition, ref. 8. Differences could result from the larger concentrations of 
H 2 /Br 2 reactants or the differences in measurement techniques - acid/base titrations for HBr 
versus laser techniques for the H/O system. Or something may be incomplete in this 
thermodynamic theory. More comparisons with data are required to completely verify a 
universal D of 7.0xl0 7 . 

One advantage of equations (11) and (12), is that the kinetic constant is continuous through 
Kp =1; and k is defined over the complete range of conditions. 
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Alternative Expression of Thermodynamic Method 

Another way of stating the thermodynamic method is simply: 

For Kp > 1: k,,=D or k, is a constant. (13) 

Equation (13) is well known for radical/radical recombinations reactions. In this work it has 
been assumed that D is the same constant for all reactions. How does the activation energy 
compare between the classical and thermodynamic methods. 

Comparison of Activation Energies between Methods 

For the classical method, the net activation energy can be computed from equation (4) and is 
given by: 

E = -R aink = nRT +E . (14) 

net 3/1 rr\ x ' 


For the thermodynamic method, E^ is computed by: 

jjn(Kp) _ AH° 

ai/T r 


50 For K p <l: E^Efy r l)RT+AH° 
For K>1: E^v, r l)RT . 


A comparison of the net activation energies at 1500 K is shown in Table 2. 


Table 2. Net Activation Energies (calories/mole) for Classical and 
Thermodynamic Methods. 


Evaluated at 1500 K 


Name 

Reaction 

-1.213 

Classical 

Thermo 

CL1 

H +02 = OH + 0 

13823. 

19039. 

CL2 

0 +H2 = OH + H 

0.062 

13619. 

2980. 

CL3 

OH +H2=H20+H 

1.544 

7980. 

2980. 

CL4 

0 + H2O=2.0OH 

-1.482 

22598. 

20043. 

CL5 

0 + 0 = 02+M 

AR= 1., H2=2.9, 02 = 1.2, H20=18.5 

10.78 

-2980. 

2980. 

CL6 

H + H = H2+M 

AR=1., H2=4.0, H=26., H20=12. 

9.509 

-2980. 

2980. 

CL7 

H + 0 = OH +M 
AR= 1., H2O=5.0 

9.571 

-1788. 

2980. 

CL8 

H + OH = H20 + M 
AR= 1., H2=2.5, H20= 16.25 

11.05 

-5961. 

2980. 

CL9 

H + 02 = H02 + M 

AR= 1., H2=3.33, 02=1.33, H20=21.3 

1.914 

-2384. 

2980. 
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CLIO 

H02+ H = OH +OH 

6.444 

1411. 

2980. 

CL11 

H02+ H = H2 +02 

7.595 

695. 

2980. 

CL12 

H02+ H = H20 + O 

7.926 

1411. 

2980. 

CL13 

H02+ O = 02 +OH 

7.657 

0 . 

2980. 

CL14 

H02+ OH =H20+02 

216.6 

0 . 

2980. 

CL15 

H02+ H02=H202 + 02 

4.740 

-1699. 

2980. 

CL 16 

H2O2 + M=2.0 OH + M 
AR=0.67, 02=0.78, H2O=6.0 

-0.215 

45502. 

52009. 

CL17 

H202+H=H02+H2 

2.855 

3775. 

2980. 

CL18 

H202 + H = H20 + OH 

10.84 

3765. 

2980. 

CL 19 

H202+0=H02 + 0H 

2.917 

6408. 

2980. 

CL20 

H202 + OH = H20 + H02 

4.400 

1431. 

2980. 


What can be stated is that when K,, is positive, the classical and thermodynamic activation 
energies are low, except for reaction CL2. For an activation energy of 2980, the kinetic 
constant increases by a factor of 2.7 when the temperature is increased from 1000 to 3000 K. 
In kinetics, this is not a large change and is represented by a classical activation energy of zero 
with an “n" of one. 

The third body coefficient for M (unity) is not counted in the k or n expression, because it was 
not included in either the free energy expression or in the IQ/Kp relation. 

Selection of Thermodynamic Reactions 

The values of the kinetic constants at 1500 K are given in Table 3 below. Individual rate 
constants for both the thermodynamic and classical methods are given. The rates included in 
the thermodynamic set, T1-T26, are shown in Table 3. Essentially all reactions which were 
possible have been written down. More reactions were included in the thermodynamic 
model than in the classical set, When K„ > 1 or log 10 K„ > 0 , the expression for the rate 
constant follows equation (12) and does not include K,, 


Table 3. Thermodynamic Reactions for the H/O System 

ORDER CLASS THERMO Evaluated at T=1500 K 


V/1VL/U 


FLF AJIAA. 

ATOM REACTIONS 


log,»K, 


log^^ 

A1 

CL6R 

T1 

M + H2 =2. OH 


-9.520 

-1.67 

0.03 

A2 

CL5R 

T2 

M + 02 =2.00 
OH REACTIONS 


-10.79 

-2.94 

-2.05 

B1 

CL7 

T3 

H + O = OH + M 


9.577 

12.94 

14.89 

B2 

CL1 

T4 

H + 02 = OH + O 


-1.213 

11.72 

11.84 

B3 

CL2 

T5 

H2 + O = OH + H 


0.062 

12.94 

12.29 

B4 



H2 + 02 =2.0OH 
H2Q REACTIONS 


-1.151 

11.79 


Cl 



2. OH + O = H20 + 

M 

20.65 

18.03 


C2 


T6 

2.0H + 02 = H20 + 

O 

9.852 

18.03 

— 

C3 


T7 

H2 + O = H20 + M 


11.12 

12.94 
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C4 

CL8 

T8 

H + OH = H20 + M 

11.06 

12.94 

15.57 

C5 

CL3 

T9 

H2 + OH = H20 + H 

1.544 

12.94 

12.66 

C6 

CL4 

T10 

0 + H20 = 2.0OH 

-1.482 

11.45 

11.11 

C7 


Til 

02 + H20 = 2.0OH + O 

-12.25 

0.711 

— 

C8 



2.0H2 + 02 = 2.0H2O 

11.46 

18.03 

— 




HQ2 REACTIONS 




D1 



H +2.00 = H02 + M 

13.09 

18.03 

— 

D2 

CL9 

T12 

H +02 = H02 + M 

1.919 

12.94 

15.31 

D3 


T13 

H2 + 2.00 = H02 + H 

3.195 

18.03 

— 

D4 

CL11R 

T14 

H2 + 02 = H02 + H 

-7.595 

5.34 

5.70 

D5 


T15 

H2+ 2.002 = 2.0HO2 

-5.66 

12.35 

— 

D6 

CLIO 

T16 

H +H02 = 2.0OH 

6.444 

12.94 

14.14 

D7 


T17 

H2 + H02 =2.0OH + H 

-2.690 

10.25 


D8 


T18 

OH + O = H02 + M 

3.492 

12.94 

— 

D9 

CL13R 

T19 

OH + 02 = H02 + 0 

-7.657 

5.28 

21.29 

DIO 

CL12 


H + H02 = H20 + O 

7.926 

12.94 

12.49 

Dll 


T20 

H2 + H02 = H20 + OH 

7.988 

12.94 

— 

D12 

CL14 

T21 

OH + H02 = H20 + 02 

9.139 

12.94 

13.30 

D13 



OH + H02 = H20 + 2.00 

-1.650 

11.28 

— 




H2Q2 REACTIONS 




El 



2. OH + 2.00 = H202 + M 

19.37 

23.11 

— 

E2 


T22 

2.0H + 02 = H202 + M 

8.578 

18.03 

— 

E3 



H2 + 2.00 = H202 + M 

9.854 

18.03 

— 

E4 



H2 + 02 = H202 + M 

-0.936 

12.00 

— 

E5 

CL 16 

T23 

M + H202 = 2.0OH 

-0.215 

7.63 

10.45 

E6 

CL17R 

T24 

H2 + H02 = H202 + H 

-2.855 

11.08 

11.68 

E7 

CL15 


2.0HO2 = H202 + 02 

4.740 

12.94 

11.27 

E8 



2.0HO2 = H202 + 2.00 

-6.050 

6.89 

— 

E9 


T25 

O + H20 = H202 + M 

-1.267 

11.67 

— 

E10 



02 + H20 = H202 +0 

-12.05 

0.878 

— 

Ell 



H + H02 = H202 + M 

6.660 

12.94 


E12 



2.0H2O = H202 + 2.0H 

-21.90 

- 8.97 


E13 

CL18 

T26 

H202 + H = H20 +0H 

10.84 

12.94 

12.48 

E14 

CL19 


H202 + O = H02 + OH 

2.917 

12.94 

12.51 

E15 

CL20 


H202 +OH = H20+ H02 

4.400 

12.94 

12.64 


TOTAL 42 REACTIONS 


The important rates were selected to complete the reaction set. The rules for rate selection 
were: 


1. The molecular/molecular reactions (B4, C8, and E4) were eliminated for two 
reasons; first, they were too fast and required an overall D value of less than 10 4 to 
obtain agreement with the data. Their reaction rate was fast not only because of a large 
kinetic constant (large K,), but also because of the high concentration of the reactants. 
Secondly, they are not included in the classical mechanism, because their rates are 
difficult to measure accurately and inclusion in the mechanism represents a parallel 
path with the free radicals. So inclusion would require a change in free radical rate 
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constants for agreement with the data. However the H0 2 molecular reactions D4 and 
D5 were included to obtain the correct ignition delay at lower temperatures. For D4 
and D5, the K,, levels are small and the forward rates are consistent with the reaction 
set. 


2. The other rates such as Cl, Dl, DIO, D13, etc., were eliminated because they were 
very small making their contribution to the overall reaction rate negligible. For 
example the combination of three radical concentrations combine to a very small 
overall rate. The subroutine ORDER and RXNTST in the LSENS program were used 
to organize the process. Subroutine ORDER orders the rates for each species at each 
output time step. RXNTST verifies that the rates are balanced and that there is no 
duplication of rates in the system. 

There are differences in several of the kinetics constants between the thermodynamic model 
and the classical model; but these differences balance out for the net system. For example, 
reaction 1 which is the initial dissociation of H2 is 50 times smaller for the thermodynamic 
method than the classical.. This appears to be not critical, because after some O is produced 
by (T2), the OH producing reaction B3 (CL2) generates H. 

Thermodynamic Kinetics Constants 


Most investigators use kinetics constants in the form of equation (4). Using the 
thermodynamic method the rate constants for all of the reactions could be computed (see 
Table 4). The values of K, were fitted to a logarithmic expression over the temperature range 
from 1000 to 3000 K. The least squares values were combined with equations (11) and (12) to 
obtain the values of A t , n t , and E l . 


Table 4. Thermodynamic Constants for the 
ORDER CLASS THERMO 


ATOM REACTIONS 


At 

CL6R 

T1 

M + H2 =2. OH 

A2 

CL5R 

T2 

M + 02 =2.00 




OH REACTIONS 

B1 

CL7 

T3 

H+ O = OH + M 

B2 

CL1 

T4 

H+ 02 = OH + O 

B3 

CL2 

T5 

H2 + 0 = OH + H 

B4 



H2 + 02 =2.0OH 




H2Q REACTIONS 

Cl 



2.0H + 0 = H20 + 

C2 


T6 

2.0H + O2 = H20 + 

C3 


T7 

H2 + O = H20 + M 

C4 

CL8 

T8 

H + OH = H20 + M 

C5 

CL3 

T9 

H2 + OH = H20 + H 

C6 

CL4 

T10 

O + H20 = 2.0OH 


H/O System 

Evaluated at T=1000 to 3000 K 


A» 

“t 

E/R 

2.71E+10 

0.982 

52575. 

6.88E+12 

0.530 

60326. 

5.74E+09 

1.0 

0 . 

2.04E+12 

0.606 

8673. 

5.74E409 

1.0 

0 . 

2.84E+12 

0.664 

9594. 

4.71E+11 

2.0 

0.0 

4.71E+11 

2.0 

0.0 

5.74E+09 

1.0 

0 . 

5.74E+09 

1.0 

0 . 

5.74 E+09 

1.0 

0.0 

8.18E+11 

0.680 

9067. 
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C7 


Til 

02+ H20 = 2.0OH + 0 

6.02E+16 

1.255 

69266. 

C8 



2.0H2 + 02 = 2.0H2O 

HQ2 REACTIONS 

4.71E+11 

2.0 

0.0 

D1 



H +2.00 = H02 + M 

4.71E+11 

2.0 

0.0 

D2 

CL9 

T12 

H +02 = H02 + M 

5.74E+09 

1.0 

0.0 

D3 


T13 

H2 + 2.00 = H02 + H 

4.71E+11 

2.0 

0.0 

D4 

CL11R 

T14 

H2 + 02 = H02 + H 

2.66E+08 

1.515 

27281. 

D5 


T15 

H2+2.0O2 = 2.0HO2 

2.62E+06 

2.05 

1986. 

D6 

CLIO 

T16 

H +H02 = 2.0OH 

5.74 E+09 

1.0 

0.0 

D7 


T17 

H2 + H02 =2.0OH + H 

3.94E+17 

-1.35 

10562. 

D8 


T18 

OH + O = H02 + M 

5.74E+09 

1.0 

0.0 

D9 

CL13R 

T19 

OH + 02 = H02 + 0 

1.91E+08 

1.457 

26359. 

DIO 

CL12 


H + H02 = H20 + 0 

5.74E+09 

1.0 

0.0 

Dll 


T20 

H2 + H02 = H20 + 0H 

5.74E+09 

1.0 

0.0 

D12 

CL14 

T21 

OH + H02 = H20 + 02 

5.74 E+09 

1.0 

0.0 

D13 



OH + H02 = H20 + 2.00 

H2Q2 REACTIQKS 

1.19E+14 

1.391 

24900. 

El 



2.0H + 2.00 = H202 + M 

3.87E+13 

3.0 

0.0 

E2 


T22 

2.0H + O2 = H202 + M 

4.71E+11 

2.0 

0.0 

E3 



H2 + 2.00 = H202 + M 

4.71E+11 

2.0 

0.0 

E4 



H2 + 02 = H202+M 

2.51E+03 

1.132 

-17282. 

E5 

CL 16 

T23 

M + H202 = 2.0OH 

7.92E+16 

-0.468 

26877. 

E6 

CL17R 

T24 

H2 + H02 = H202 + H 

2.09E+07 

1.600 

8012. 

E7 

CL15 


2.0HO2 = H202 + 02 

5.74E+09 

1.0 

0.0 

E8 



2.0HO2 = H202 + 2.00 

4.44E+13 

1.614 

41057. 

E9 


T25 

0 + 1120 = H202 + M 

7.23E+02 

1.15 

-17810. 

E10 



02 + H20 = H202 +0 

7.11E+07 

1.68 

42516. 

Ell 



H + H02 = H202 + M 

5.74E+09 

1.0 

0.0 

E12 



2.0H2O = H202 + 2.0H 

7.93E+09 

2.68 

94564. 

E13 

CL18 

T26 

H202 + H = H20 +OH 

5.74E+09 

1.0 

0.0 

E14 

CL19 


H202 + 0 = H02 + OH 

5.74E+09 

1.0 

0.0 

E15 

CL20 


H202 +OH = H20+ H02 

5.74E+09 

1.0 

0.0 


TOTAL 42 REACTIONS 


When Kp is greater than one, many rates are computed as the same values. However the net 
forward rate depends on the concentrations as well as the kinetic constants. 

In the next section the results from the classical method and from the thermodynamic method 
using equations (11) and (12) to compute the kinetics for all reactions are compared with 
experimental data. 
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Comparison of Both Methods to Experimental Data 

The LSENS code was modified to compute the forward kinetic constant using the Gibbs free 
energies. One could also use the constants in Table 4. The code already had provisions to 
calculate the reverse rate from the equilibrium constant, equation (1). Only about five lines 
had to be added to the code in order to compute the forward rate from the thermodynamic data 
and equations (11) or (12). 

The thermodynamic method was first tested using the classical mechanism equations (CL1- 
CL20 but with new constants). The experimental data is given in terms of the OH absorption, 
which is related to concentration by: 

AB OH =l-exp(1.5xl0 7 L C OH moIes/cc); 
where L=path length in cm. 

The agreement between the experimental data and the classical mechanism was excellent for 
the single test case shown in Figure 1 . This single test case was used for the thermodynamic 
method with the constant D of 3x10* as reported in reference 1. The classical values of M 
were used. The results show the predicted ignition time was four times too fast. 



Figure 1: Comparison of Classical and Thermodynamic Predictions to Experimental data, 
T=1556 K, P= 0.751 atm, XH2=0.02, X02=0.002, or Phi =2. 


When D was lowered to 7xl0 7 with the classical values of M , the agreement was much better. 
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However, the thermodynamic method did not provide any method for 
computing M. 

The most confusing thing about third body efficiencies is the variety of values given in the 
literature. For instance, for the reaction H 2 + M = 2 H + M, reference 2 give values of M 
as H 2 =4, H=26. , and H 2 0= 12. While reference 9 reports M values of H 2 =0. 14T 04 , 

H 2 0=93.8 T 025 , C 0=1.6, N 2 =1.6, 0 2 =1.6, and C 0 2 =858/T. And reference 10 gives 
values of H 2 =4.1, H 2 0=15., 0 2 =2., and N 2 =2. In this report, in order to keep things as 
simple as possible, all third body efficiencies M were set to unity and all collisions assist in 
the reaction - even self collisions. So the concentration of the third body M is given by: 
C m = P ,m w- 

However to check this approach, computations were made without M; but the fit over the 
complete temperature and pressure range of the data was poor. It is logical to include M in 
the set in order to get the correct pressure and temperature dependence for monomolecular 
reactions. All of the efficiency M values were taken as unity. This is consistent with the 
free energy formulation from equation (6) and (7), where the concentration of the diluent 
would affect the mixing term £ nj (In p ; ) and would not cancel out when there was a change 
in the number of moles for the reaction. 

There is excellent agreement between the shapes of the OH rise curves for the classical and 
thermodynamic methods. Good agreement of the OH overshoot was obtained between the 
classical and thermodynamic methods 

With the thermodynamic approach all of the kinetic constants are calculated by the LSENS 
program using thermodynamic data, so one only needs to write down the chemical reaction to 
include it in the system. 

The final comparisons are shown in Table 5. The times t^, and t^, are the reaction times 
required to reach 25, 50, and 75 percent of the maximum OH concentration. This table 
represents the minimum and maximum of the temperatures, and pressures for the 9 different 
compositions run in ref 2. A total of 200 cases were reported in ref. 2 but only 20 are selected 
here. 
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Table 5. Comparison of Classical and Thermodynamic Methods to 
Experimental Data. *S= Single path length, L = 6.35 cm. 

D= Double path length, L=12.7 cm 





7D 

1525. 

8D 

2409. 

9D 

1527. 

10D 

2211. 

ns 

1243. 

12D 

2414. 

13S 

1509. 



3.770 0.002 0.0005 Exper 573. 

<{>=2 Class 635. 
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Experimental times t50- microseconds 


Figure 2: Predictions versus Experimental 

The classical predictions are the solid square symbols and the thermodynamic 
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predictions are the open triangle symbols in Figure 2 for the t^ times from 
Table 5 above. The classical predictions agree with the experiments to ± 10 % while the 
thermodynamic predictions agree with the data to ±30%. This is good agreement for the 
thermodynamic method considering it contains only a single constant D. 

Reactions CL1, CL2, and CL3 are the most important reactions of the system. CL1 and CL2 
produce OH and CL3 consumes OH to produce H 2 0. The ratio of the thermodynamic to 
classical predictions for these three kinetic constants are shown in Figure 3. Perfect 
agreement between the two methods would result in a ratio of one. 



Temperature K 


Figure 3: Comparison of Thermodynamic and Classical Rates, D=7.0 E 7. 


The discrepancy is greater at the lower temperature than at a higher temperature. This is 
because the free radical concentrations are low at low temperatures, making measurement 
accuracy low. Also at low temperature the rates are slow so that systematic errors are 
magnified on a ratio scale. In addition if an error is present in one rate, it can propagate to 
other reactions because of the coupling effect. So the higher consumption rate of CL2 is 
compensated by the low rate of CL1 and the high consumption rate of CL3. 
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Reevaluation of the Hydrogen-Bromine System 


From the above study and the inclusion of the third body efficiencies in the rate equations, 
the hydrogen-bromine system from ref. 1 was recomputed to determine if there was still a 
difference in the value of D between the two systems. The rate equations used for this study 
were: 


1. M + Br 2 = 2 Br + M 

The equations in 

T1 

Br 2 =2 Br 

2. M + H 2 = 2H + M 

ref. 1 were: 

T2 

H 2 = 2 H 

3. Br + H 2 = HBr + H 


T3 

0.5H 2 +Br=HBr 

4. H + Br 2 = HBr + Br 2 


T4 

H+0.5Br 2 =HBr 

5. M + H + Br = HBr -1- M 


T5 

H+Br= HBr 


Reference 1 treated the “mechanism” equations as simply mass balance equations between 
species. This report is treating the “mechanism” as all possible interactions between species 
with the path determined by the Gibbs free energy gradient. The new equations are essentially 
the classical mechanism or the thermodynamic equations without the previously used 
fractional stoichiometric coefficients. The molecule/molecule reaction between H 2 and Br 2 
(H 2 +Br 2 =2 HBr) was not included. The best value of D for comparison to the HBr data was 
still 3xl0 8 , which is the same as determined in ref.l. The values of the new predicted results 
are shown in Table 6. 


Table 6. HBr Mole Fraction Data at 574.3 K ; p=1.8 atm. 
Initial Conditions: Xh 2=0.657 X Brl =0.343 


Time 

Data 

Classical 

Thermodynamic 

seconds 

^HBr 

^HBr 

^HBr 

870 

0.1559 

0.1404 

0.1855 

1470 

0.2295 

0.2165 

0.2580 

2070 

0.2940 

0.2797 

0.3106 

2970 

0.3830 

0.3559 

0.3668 

4770 

0.4877 

0.4613 

0.4389 

5970 

0.5373 

! 0.5086 

0.4710 

7470 

0.5829 

0.5512 

0.5015 

8970 

0.6102 

0.5814 

0.5240 

Eauilibrium 

0.686 
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The thermodynamic method is initially 19% too high, crosses the experimental data 
at a time of 2440 seconds and finishes 14% too low. Decreasing the value of D would 
decrease the concentration for all times and shift the crossover point to an earlier time. A value 
of D of 3x10® was considered to give good agreement with the data. The classical values are 
within 10% for the whole time series. 

Comparison of Other Hydrogen-Bromine Data with Model 

Now that D is fixed, the computed results can be compared with other data. In general the 
prediction of the effects of temperature was much better than the results of ref. 1. These 
calculations included the effect of third bodies M, and did not use fractional coefficients in the 
reaction set. The calculations of the effect of increased HBr concentration, reduced 
temperatures, and changing initial conditions compared well. It is possible that the difference 
in the value of D for the two systems could be a result of the measurement techniques, 
impurities, or the large concentrations of the reactants. The HBr data was obtained using 
acid/base titration techniques with reaction times of minutes compared to the H/O system 
which used laser absorption with reaction times of microseconds. As yet unpublished work on 
the CjHy/CO/NO system suggests that in general 7xl0 7 works well for D. 

Conclusions 

The thermodynamic method adequately predicted the reaction rate constants for the 
Hydrogen/Oxygen system. The thermodynamic method is based on the gradient of the Gibbs 
free energy and a single proportionality constant D to determine the kinetic rate constants. 
Using this method the rate constants for any gas phase reaction can be computed from 
thermodynamic properties. Any gas phase kinetic constant can be computed given the 
thermodynamic data from the equations: 

k cf =DK p (RT) Sv,rl for K p < 1 (11) 

or k c =D(RT) SVfi 1 for K p > 1. (12) 

The deviation between predictions and experimental data for the thermodynamic method was 
greater than for the classical method, but there is only one constant D (7xl0 7 )in the 
thermodynamic method to determine the whole H/O system. All of the third body efficiencies M 
are taken to be unity. This new thermodynamic based kinetics model appears very simple; but 
it appears to work. There was a larger deviation with experimental data, however there is 
probably error in the experimental data which cannot be compensated for with only one 
constant. One reason that this method has been shown to work is that the experimental data is 
very good, so that this thermodynamic method is able to match with it. 
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